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ABSTRACT 


Lack  of  measurements  has  led  to  a  poor  understanding  of  tidal  eurrents  within  the  Mouth 
of  the  Columbia  River  (MCR).  Shipboard  ADCP  measurements  are  a  teehnically  viable 
means  to  measure  these  eurrents  and  their  vertical  structures.  Shipboard  ADCP 
measurements  from  a  three-week  long  eruise  conducted  by  the  RA^  Point  Sur  are  used  in 
this  study.  Transects  along  the  main  shipping  channel  and  aeross  the  bar  were  extraeted 
using  eriteria  of  eonstant  speed  and  heading.  Vertical  sections  from  each  transect  were 
compared  with  predictions  from  the  SELFE  model  in  order  to  evaluate  the  models’ 
ability  to  predict  these  currents. 

Observed  flood  currents  within  the  river  are  bottom  intensified  with  slight  river 
outflow  at  the  surfaee  early  in  the  flood  cycle.  Complete  reversal  of  the  flow  throughout 
the  column  is  experienced  later  as  the  flood  fully  develops.  Ebb  currents  are  mueh 
stronger  and  more  homogenous  throughout  the  water  eolumn  due  to  the  combination  of 
river  outflow  and  tidal  current. 

Inside  the  river  mouth,  the  SEEFE  model  agrees  qualitatively  well  with 
observations,  with  only  a  few  diserepancies.  Outside  the  river  mouth,  the  SEEFE  model 
does  not  resolve  the  observed  vertical  structure,  possibly  due  to  coarse  grid  resolution  and 
limited  observational  input. 
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I.  INTRODUCTION 


A.  MOTIVATION  AND  BACKGROUND 

Coastal  environments  are  vital  areas  of  interest  to  all  aspeets  of  soeiety.  Important 
aetivities  in  the  eoastal  zone  inelude  import  and  export  of  goods  as  well  as  eommereial 
fishing,  reereation  on  beaehes  and  waterways,  U.S.  Coast  Guard  seareh  and  reseue  and 
naval  maritime  operations,  energy  produetion  and  eonstruetion  of  eoastal  struetures 
(piers,  jetties,  and  bars).  The  study  of  these  environments  is  important  beeause  a  better 
understanding  of  their  dynamies  leads  to  safer  navigation,  more  aeeurate  ehart 
produetion,  safer  design  and  building  of  eoastal  struetures,  and  a  better  understanding  of 
the  evolution  of  eoastal  regions  through  eurrents  and  sediment  transport. 

One  important  aspeet  of  eoastal  researeh  is  the  study  of  tidal  eurrents.  The 
strength  and  orientation  of  these  eurrents  affeet  other  important  aspeets  of  eoastal 
environments  sueh  as  sediment  transport,  bathymetry,  and  wave  interaetions  at  river 
mouths.  The  ability  to  aoeurately  prediet  the  behavior  of  eurrents  throughout  a  tidal  eyele 
is  a  key  step  toward  safety  of  navigation,  and  the  management  of  eoastal  regions. 

This  study  foeuses  on  the  eurrents  in  the  Mouth  of  the  Columbia  River  (MCR) 
(Figure  1).  The  Columbia  River  is  the  largest  river  on  the  North  Ameriean  Paeifie  Coast. 
The  river’s  strong  outflow,  eoupled  with  its  strong  tidal  eurrents,  ereate  an  estuary  at  the 
mouth  of  the  river  with  very  dynamie  oeeanographie  features,  eharaeterized  by  primarily 
broad  shoals  separating  narrowing  ehannels.  Atmospherie  proeesses  are  seeondary  in 
nature  to  the  effeets  of  river  flow  and  tidal  infiuenee.  Tidal  eurrents  in  partieular,  are  a 
dominant  faetor  in  determining  the  flow  phenomena  oeeurring  within  the  estuary  (Jay  and 
Smith  1990). 
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Figure  1 .  Mouth  of  the  Columbia  River  with  the  North  and  South  Jetties 

labeled,  (from  Wikipedia  cited  2014). 


1,  Tides 

Tides  are  the  periodic  rise  and  fall  of  the  sea  surface  caused  by  gravitational 
attractive  forces  between  the  sun  and  moon  and  the  rotating  Earth.  The  moon’s 
gravitational  pull  is  the  principle  tide-generating  force  due  to  the  sun’s  greater  distance 
from  Earth  (Bowditch  2002). 

Coastal  regions  all  over  the  world  experience  one  of  three  different  types  of  tides; 
diurnal,  semi-diurnal,  and  mixed  tides  (Gross  2005).  Diurnal  tides  occur  when  a  coastal 
region  experiences  only  one  high  tide  and  one  low  tide  per  day  (Gross  2005). 
Semidiurnal  tides  occur  when  a  coastal  region  experiences  two  high  tides  and  two  low 
tides  per  day.  Semi-diurnal  tides  are  observed  at  the  equator  year  round  (Gross  2005). 
Mixed  tides  occur  when  a  coastal  region  experiences  two  unequal  high  tides  and  low 
tides  within  a  day. 
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The  tidal  range  is  affected  by  the  alignment  of  the  Earth-sun-moon  system 
(Bowditch  2002).  When  the  Earth,  sun  and  moon  are  aligned  the  lunar  and  solar  tractive 
forces  combine  to  produce  the  greatest  tidal  ranges  of  the  month  called  spring  tides 
(Eigure  2A).  The  sun  and  moon  at  right  angles  to  the  Earth  result  in  the  lowest  tidal  range 
of  the  month,  called  neap  tides  (Eigure  2B.) 


Eigure  2.  (A)  Spring  tides  occur  when  the  solar  and  lunar  tractive  forces  are 

aligned  with  the  Earth.  (B)  Neap  tides  occur  when  the  solar  and  lunar 
tractive  forces  are  at  right  angles  to  each  other  (from  Bowditch  2002). 


2,  Tidal  Currents  in  Coastal  Regions 

In  coastal  regions  where  flow  is  restricted  by  bathymetry  such  as  channels  or 
straits,  tidal  currents  alternate  between  flood  and  ebb,  filling  and  draining  the  bay  or 
estuary  (Eigure  3).  During  slack  water,  the  flow  reversal  between  ebb  and  flood,  the  tide 
induced  flows  are  very  weak. 
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Rotary  tkJal 
currents 


Figures.  Flood  and  ebb  currents  within  a  harbor.  Flood  current  flows 
upstream  while  ebb  current  flows  downstream  (from  Gross  2005). 


Like  the  tidal  sea  level  variations,  tidal  currents  are  periodic  in  nature,  with 
variations  in  amplitude  over  the  tidal  cycle.  The  strongest  tidal  currents  occur  during 
spring  tides  while  the  weakest  tidal  currents  occur  during  neap  tides  (Gross  2005).  This 
periodicity  exhibited  by  tidal  currents  distinguishes  them  from  other  non-tidal  currents. 
Other  variations  in  tides  and  tidal  currents  occur  as  a  result  of  the  moon’s  continuous 
changes  in  phase  (29.5  days),  distance  from  the  Earth  (27.5  days),  and  declination  to  the 
Earth  (27.3  days)  (Disney  and  Overshiner  1925).  Changes  in  the  moon’s  phase  and 
distance  from  Earth  cause  changes  in  tidal  current  strength  that  are  approximately 
proportional  to  changes  in  tidal  range.  Changes  in  the  moon’s  declination,  however,  lead 
to  changes  in  tidal  current  strength  that  are  approximately  half  of  the  changes  in  tidal 
range  (Disney  and  Overshiner,  1925). 
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3. 


Non-Tidal  Currents 


Non-tidal  constituents  in  the  form  of  river  outflow  and  wind-driven  eurrents 
usually  aeeompany  tidal  eurrents  in  eoastal  regions.  Wind  drives  surfaee  eurrents  that 
deerease  with  depth  due  to  seabed  frietion.  Temperature  differenees  between  two 
separate  bodies  of  water,  as  well  as  density  differenees  between  freshwater  and  salt 
water,  eause  horizontal  stratifieation  within  the  water  eolumn  that  leads  to  density  driven 
currents  (Koeh  and  Sun  1999). 

4,  Time  Scales  of  Tidal  and  Non-Tidal  Currents 

Tidal  and  non-tidal  current  time  scales  are  shown  in  Figure  2.  Density  driven 
currents  operate  over  many  different  time  scales  depending  on  the  forcing  mechanisms. 
Ebb  tides  can  transport  fresh  water  into  estuaries  or  out  onto  continental  shelves  over 
short  time  periods.  Flash  floods  resulting  from  rainfall  can  last  for  days.  Seasonal 
processes  such  as  runoff  and  evaporation  can  last  as  long  as  the  season  in  which  they  are 
happening,  weeks  or  more.  Tidal  and  wind  driven  currents  in  coastal  areas  often  have 
similar  time  scales  because  of  the  diurnal  cycle  of  the  sea  breeze. 
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5,  Mouth  of  the  Columbia  River 

a.  The  Columbia  River 

The  Columbia  River  drainage  basin  eovers  an  area  approximately  660,480  km^  in 
seven  U.S.  states  and  two  Canadian  provinces  (Simenstad  et  al.  2002).  With  an  average 
discharge  of  approximately  7,300  m^sec’',  the  Columbia  River  is  the  largest  river  on  the 
U.S.  Pacific  coast  and  the  fourth  largest  river  in  the  U.S.,  by  volume  (Jay  et  al.  1984). 
The  Columbia  River  originates  in  the  Rocky  Mountains  of  British  Columbia  and  flows  in 
a  mean  southern  direction  to  empty  out  into  the  Pacific  Ocean  along  the  Washington- 
Oregon  border  in  the  northwest  U.S.  The  river’s  large  flow  rate,  combined  with  a  steep 
gradient  from  its  mountainous,  decent  gives  the  Columbia  River  tremendous  potential  for 
generating  hydroelectric  power.  The  river  has  been  dammed  hundreds  of  times  in  the  last 
century  for  hydroelectric  power  generation,  irrigation,  and  flood  control  (Elias  et  al. 
2012). 
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b.  Bathymetric  Features  of  MCR  and  Offshore  Area 

The  Columbia  River  bar  extends  offshore  from  the  mouth  of  the  river 
approximately  8  km.  It  is  a  relatively  flat,  and  shallow  (~20  m)  region  that  is  radially 
symmetric  about  the  tip  of  the  North  Jetty  (Kilcher  and  Nash  2010).  The  estuary  is 
constricted  at  the  mouth  of  the  river  by  two  jetties:  the  North  Jetty  located  at  the  tip  of 
Peacock  Spit  and  the  South  Jetty  located  at  the  tip  of  Clatsop  Spit.  The  river  narrows 
from  10  km  within  the  estuary  to  3.5  km  at  the  mouth  of  the  river  before  flowing  into  the 
Pacific  Ocean  (Kilcher  and  Nash  2010).  The  estuary  is  also  rather  shallow  except  for  the 
two  shipping  channels  that  cut  through  the  bar  and  run  through  the  estuary.  The  North 
Channel  splits  from  the  main  channel  and  touches  on  the  southern  coast  of  Washington 
within  the  estuary.  The  main  Navigation  Channel  hugs  the  southern  edge  of  the  estuary 
and  continues  upriver  southeast.  Both  shipping  channels  average  20  m  depth  and  are  the 
deepest  points  within  the  estuary  (Figure  5). 
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c.  Tides  and  Tidal  Currents  within  the  MCR 

Tidal  forcing  dominates  the  dynamics  within  the  estuary,  but  river  flow  takes  on 
an  increasingly  important  role  upriver  from  the  entrance.  The  Columbia  River  is  affected 
by  tides,  as  far  as  243  km  upriver  from  the  mouth  (Templeton  and  Jay  2013).  Tides 
within  the  MCR  are  mixed,  consisting  of  two  high  waters  and  two  low  waters  over  a  24.8 
hour  period,  with  a  mean  diurnal  range  of  approximately  1.8  m  (Jay  et  al.  2011).  At  the 
mouth  of  the  river  the  mean  and  diurnal  tidal  ranges  are  2  m  and  2.6  m,  respeetively 
(Lutz  et  al.  1975).  Figure  6  depiets  NOAA  predieted  tides  at  the  Hammond,  OR  station 
for  the  period  of  25  May  2013  to  01  Jun  2013.  This  is  a  typical  mixed,  semi-diurnal  tidal 
range  experieneed  in  this  eoastal  region. 
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Figure  6.  NOAA  predieted  tides  for  the  period  of  25May2013-01  Jun2013 

(from  NOAA  Tides  cited  2014). 


Tidal  currents  are  stronger  within  the  North  Channel  because  most  of  the  tidal 
prism  of  the  lower  estuary  is  filled  from  the  North  Channel  (Jay  and  Smith  1990).  Much  of 
the  tidal  transport  is  applied  to  fill  the  tidal  prism  over  the  estuary,  and  therefore,  tidal 
eurrent  veloeities  deerease  strongly  upriver  (Jay  and  Smith  1990).  During  the  ebb  cycle, 
surfaee  eurrents  ean  routinely  reaeh  up  to  and  exceed  2-3  ms"'  causing  a  significant  increase 
in  the  height  and  steepness  of  opposing  surfaee  waves  in  a  matter  of  less  than  2-3  hours 
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(Gonzales  2000).  In  regions  where  the  water  eolumn  is  well  mixed,  bottom  friction  slows 
the  currents  greatly,  leading  to  velocities  decreasing  with  depth  from  the  surface. 

Strength  of  the  river  flow  also  tends  to  have  an  effect  on  tidal  currents.  When 
non-tidal  constituents  such  as  river  flow  and  wind  forcing  are  not  considered,  ebb  and 
flood  periods  and  magnitudes  are  approximately  equal.  The  strength  of  the  Columbia 
River  flow,  however,  acts  to  increase  the  strength  of  the  ebb  current  while  decreasing  the 
strength  of  the  flood  current.  River  flow  also  impacts  the  length  and  time  period,  of  the 
tidal  cycles,  causing  ebb  duration  to  increase  and  flood  duration  to  decrease  (Bowditch 
2002). 

B.  NAVAL  RELEVANCE 

The  Columbia  River  bar  is  one  of  the  major  landmarks  of  the  infamous 
“Graveyard  of  the  Pacific”.  This  region  is  characterized  by  intense  wave-current 
interactions,  and  has  been  recognized  as  one  of  the  most  dynamic  and  dangerous  coastal 
regions  in  the  entire  world.  With  an  annual  passage  of  approximately  2,000  vessels, 
safety  of  navigation  through  the  channel  has  become  paramount  (Simenstad  et  al.  2002). 

Cape  Disappointment,  located  just  north  of  Peacock  Spit,  is  home  to  Coast  Guard 
Station  Cape  Disappointment.  Gonzales  (1984)  notes  that  statistics  from  this  station 
showed  that  an  average  of  850  search-and-rescue  missions  were  conducted  each  year. 
The  construction  of  the  North  and  South  jetties  at  the  entrance  of  the  river  greatly 
improved  navigation  but  shifted  the  morphology  of  the  inlet.  Because  of  this,  dredging  is 
required  to  maintain  the  channel  depth  at  its  required  17  m.  The  jetties  require  frequent 
maintenance  due  to  severe  wave  and  tide  interactions  (Elias  et  al.  2012). 

Safety  of  navigation  is  important  through  this  region  because  of  the  dynamic 
environmental  conditions,  but  also  because  of  the  regional  economic  situation.  The  ports 
of  Portland  and  Vancouver  are  accessed  via  the  Columbia  River.  Research  from  Jay  et 
al.  (2011)  reveals  that  ongoing  changes  in  water  levels  and  tidal  ranges  affect  ships  draft 
within  the  channel  as  well  as  shallow-water  habitat.  Decreasing  water  levels  coupled 
with  increasing  ship  drafts  lead  to  an  increase  in  difficulty,  and  risk  navigating  the 
Columbia  River. 
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C.  SCOPE 

The  goals  of  this  research  are  twofold;  characterize  the  vertical  structure  of  tidal 
currents  within  the  MCR  over  various  tidal  cycles,  and  validate  the  SELFE  model’s 
ability  to  predict  these  currents. 

Characterizing  the  vertical  structure  of  the  tidal  currents  is  important  because  of 
the  strong  stratification  typically  found  in  river  mouths.  Saline  water  from  the  ocean  is 
denser  than  fresh  water  from  upriver,  so  this  difference  in  density  causes  strong  shear 
throughout  the  column  as  these  two  water  masses  meet  (Elias  et  al.  2012).  The  ADCP 
measurements  collected  onboard  the  R/V  Point  Sur  in  the  ONR  RIVET  II  Experiment 
provides  an  opportunity  to  investigate  the  vertical  structure  of  tidal  currents  in  the  MCR. 

Despite  the  significant  economic  and  social  importance  of  the  MCR,  relatively 
few  detailed  measurements  of  the  tidal  flows  have  been  reported.  As  a  result,  little  is 
known  about  its  vertical  tidal  structure. 

Evaluating  operational  models  of  the  MCR  circulation  such  as  the  SEEFE  model, 
tested  in  this  research,  will  provide  mariners  with  a  more  accurate  depiction  of  the  tidal 
current  variability  within  the  river  mouth.  Understanding  these  currents  and  their 
structure,  and  helping  to  validate  a  model  that  can  eventually  be  used  as  a  regular 
navigation  tool  is  a  step  in  the  direction  to  safer  navigation  within  the  estuary. 
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II.  MODELING  THE  COLUMBIA  RIVER 


A.  SELFE 

There  are  a  few  models  currently  in  use  to  predict  hydrodynamic  conditions 
within  the  Columbia  River  estuary.  The  Semi-implicit  Eulerian-Lagrangian  Finite 
Element  model  (SEEFE)  is  an  operational  modeling  system  that  is  a  refinement  of  the 
well-established  Eulerian-Eagrangian  Circulation  model  (EECIRC)  (Burla  et  al.  2010a). 
This  semi-implicit,  unstructured  grid  model  solves  the  3D  barotropic  shallow-water 
equations  and  can  be  used  to  predict  3D  velocities,  water  levels,  temperature,  and  salinity 
at  each  point  in  space  and  time  over  the  entire  computational  domain  (Burla  et  al.  2010b). 
SEEFE  utilizes  a  finite-element  framework  with  a  terrain-following  coordinate  system 
approaching  the  bottom  boundary  layer  (Burla  et  al.  2010b). 

The  SEEFE  model  came  about  by  efforts  to  improve  older  unstructured  grid 
models  such  as  EECIRC  (Zhang  and  Baptista  2008).  SEEFE  retains  many  important 
base  features  from  previous  models  to  include  computational  efficiency  from  semi- 
implicit  Eulerian-Eagrangian  finite  volume  modeling,  but  allows  for  superior  flexibility 
in  representing  bathymetry  by  relaxing  restrictions  on  grids,  and  using  higher-order  shape 
functions  to  calculate  elevations  (Zhang  and  Baptista  2008).  SEEFE  can  utilize  either  S 
coordinates  or  hybrid  SZ  coordinates  to  represent  a  vertical  grid  that  resembles 
unstructured  grids,  which  allow  for  superior  representation  of  bathymetry  and  vertical 
structure  of  the  water  column  (Zhang  and  Baptista  2008).  SEEFE  solves  for  3D  velocity, 
3D  temperature  and  salinity,  and  free-surface  elevation  (Zhang  and  Baptista  2008). 
Figure  7  shows  the  MCR  SEEFE  model  domain  grid. 
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Figure  7.  (a)  Entire  MCR  SELFE  model  grid,  (b)  MCR  estuary  SEEFE  model 

grid,  and  (e)  MCR  tidal  freshwater  SEEFE  model  grid  (after  Kama 

2014). 


B.  SELFE  PHYSICS 

SELFE  solves  for  the  3D  shallow-water  equations  by  using  a  eombination  of  the 
hydrostatie  and  Boussinesq  approximations,  as  well  as  transport  equations  for  salt  and 
heat  (Zhang  and  Baptista  2008).  The  Generic  Length  Scale  turbulence  closure  of  Umlauf 
and  Burchard  (2003)  is  used  to  calculate  turbulence  within  SELFE. 

Most  ocean  models  perform  poorly  at  resolving  the  bottom  boundary  layer  when 
using  the  no-slip  condition,  (m=w=0).  SELFE,  however,  uses  a  balance  between  the 
internal  Reynolds  stress  and  the  bottom  frictional  stress  to  resolve  the  bottom  boundary 
layer  (Zhang  and  Baptista  2008).  Different  boundary  layers  (i.e.,  laminar  vs.  turbulent) 
may  require  different  forms  of  bottom  stress  (Zhang  and  Baptista  2008).  SELFE  is 
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flexible  enough  to  be  applied  to  various  boundary  layer  types  and  bottom  stresses  (Zhang 
and  Baptista  2008). 

C.  SELFE  NUMERICS 

SELFE  has  several  advantages  over  its  predecessors  that  enhance  numerical 
efficiency  and  accuracy.  First,  SEEFE  solves  the  shallow  water  equations  system 
described  above  using  a  finite-element  and  finite -volume  scheme  (Zhang  and  Baptista 
2008).  This  eliminates  the  requirement  for  splitting  the  internal  and  external  modes,  thus 
eliminating  errors  that  are  related  with  mode  splitting  (Shchepetkin  and  McWilliams 
2005).  Second,  advection  terms  within  the  momentum  equations  are  treated  using  an 
Eulerian-Eagrangian  method  which  reduces  numerical  stability  constraints  (Zhang  and 
Baptista  2008).  Finally,  SEEFE  uses  a  semi-implicit  scheme  to  all  equations,  while 
solving  the  continuity  and  momentum  equations  simultaneously,  eliminating  stability 
restrictions  associated  with  previous  models  (Zhang  and  Baptista  2008). 

D.  IMPLEMENTATION  WITHIN  THE  MCR 

Within  the  MCR,  SEEFE  is  used  to  produce  data  via  the  application  of  the  model 
to  a  variety  of  inputs.  River  discharges  at  Bonneville  Dam  are  from  the  U.S.  Army  Corps 
of  Engineers,  while  other  river  inputs  are  from  U.S.  Geological  Survey  (USGS) 
instrumentation  throughout  the  river  and  estuary  (Turner  2014). 

Atmospheric  conditions,  air  temperature,  solar  flux,  and  winds  are  taken  from  the 
National  Centers  for  Environmental  Prediction  (NCEP)  North  American  Regional 
Reanalysis  (NARR).  This  database  contains  data  sets  from  numerical  weather 
predictions  to  include  model  outputs,  weather  model  outputs  and  inputs,  as  well  as 
climate  model  data  sets  (Turner  2014). 

Ocean  boundary  conditions  for  salinity,  temperature,  and  sea  surface  height  come 
from  two  different  models  depending  on  the  dates  of  the  data  (Turner  2014).  Data  from 
years  prior  to  2013  are  acquired  from  the  Naval  Research  Eaboratories  Navy  Coastal 
Ocean  Model  (NCOM).  This  model  has  since  been  retired  and  all  subsequent  data  come 
from  Hybrid  Coordinate  Ocean  Model  (HYCOM)  (Turner  2014). 


13 


Tides  are  provided  to  SELFE  via  the  Advaneed  Cireulation  model  (ADCIRC) 
(Turner  2014).  ADCIRC  is  a  finite  element  hydrodynamie  model  based  on  unstruetured 
grids  similar  to  SEEFE  (Meyers  and  Baptista  2001).  ADCIRC  applies  variable 
resolution,  whieh  allows  the  model  to  transition  between  a  eoarser  representation  in 
deeper  oeean  environments  and  a  finer  representation  in  nearshore  environments  (Meyers 
and  Baptista  2001). 

Bathymetry  within  SEEFE  is  a  eomposite  of  multiple  sources.  Bathymetry  is 
loaded  to  the  model  grid  sequentially  from  the  following  progressively  higher  resolution 
data  sets: 

•  NGDC  ET0P02v2 

•  NGDC  3arc  second  Coastal  Relief  Model 

•  NGDC  1/3  arc  second  raster  for  EaPush,  Taholah,  and  Astoria 

•  DOGAMI  combined  bathymetry  and  LiDAR  dataset  for  the  Columbia 
River  Estuary  up  to  Bonneville  and  Willamette  Falls 

•  USACE  2012  bank  to  bank  survey  (from  Jetty  A  to  Tongue  Point) 

•  USACE  yearly  cross  channel  surveys  for  the  area  around  the  navigation 
channel 

•  USACE  monthly  along  channel  surveys  for  the  navigation  channel. 

Figure  8  shows  the  MCR  SEEFE  bathymetry  grid. 
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Figure  8.  Bathymetry  for  (a)  Entire  MCR  SELFE  model  grid,  (b)  MCR  estuary 
SEEFE  model  grid,  and  (c)  MCR  tidal  freshwater  SEEFE  model  grid 

(after  Kama  2014). 
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III.  FIELD  EXPERIMENT 


A.  DATA  COLLECTION  AND  PREVIOUS  RESEARCH 

Data  for  this  research  were  collected  during  the  Office  of  Naval  Research  (ONR) 
sponsored  RIVET  II  Experiment  in  the  MCR.  The  focus  of  this  project,  involving 
scientists  from  many  institutions,  is  the  hydrodynamic  and  morphodynamic  processes 
resulting  from  waves,  tidal  currents  and  river  outflow  in  this  complex  estuary.  The  NFS 
team  in  collaboration  with  Thesis  Research  (Dr.  Tim  Janssen)  was  responsible  for 
collecting  in-situ  wave  and  current  measurements. 

The  data  used  in  this  research  were  collected  during  a  research  cruise  within  the 
MCR  from  25  May  2013  to  14  Jun  2013  aboard  the  scientific  research  vessel  lUV  Point 
Sur.  Pearman  (2014)  researched  wave  evolution  within  river  mouths  and  tidal  inlets  via 
data  collected  by  a  fleet  of  Wave  Resolving  Drifters  (WRDs)  that  were  deployed  and 
recovered  several  times  over  the  course  of  the  cruise.  Three  slightly  different  variants  of 
the  drifters  were  used  to  measure  currents  and  waves  via  various  onboard  arrangements 
of  equipment  to  include  GPS  loggers,  accelerometers,  and  Inertial  Motion  Units  (IMU). 
The  drifters  measured  horizontal  orbital  velocities  in  the  swell  band,  as  well  as  vertical 
accelerations  of  higher- frequency  wind  waves  (Pearman  2014). 

Porte  11  (2013)  validated  the  use  of  IMUs  as  an  accurate,  and  low-cost  alternative 
for  use  in  recording  surface  waves.  Data  collected  from  the  WRDs  deployed  during  this 
experiment  were  also  used  to  help  validate  the  results  Portell  (2013)  concluded  in  his 
experiments.  Portell  (2013)  constructed  a  wave  simulation  device  to  test  IMUs  at 
multiple  wave  frequencies  and  unit  configurations  within  a  laboratory  setting.  These 
results  were  compared  to  observations  obtained  from  IMUs  deployed  during  the 
experiment  to  estimate  error  analysis  and  accuracy  (Portell  2013). 

Finally,  ADCP  data  were  collected  from  both  a  shipboard  mounted  ADCP  unit,  as 
well  as  a  bottom-moored  ADCP  unit.  Data  from  these  units  were  collected  at  a  rate  of 
one  profile  every  two  minutes.  The  shipboard  ADCP  data  will  be  used  here  to 
characterize  the  vertical  structure  of  the  tidal  currents  within  the  MCR.  These  data  will 
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also  be  compared  to  the  SELFE  model  representations  of  the  MCR  in  order  to  validate 
the  accuracy  of  the  model  within  this  region. 

B,  FIELD  SITE 

All  of  the  data  used  in  this  research  were  collected  in  the  main  Navigation 
Channel  located  in  the  southern  portion  of  the  Columbia  River  estuary,  and  just  outside 
of  the  estuary  off  the  bar  (Figure  9).  Transects  inside  the  mouth  of  the  river  were 
collected  in  the  region  of  the  estuary  between  the  entrance  (North  and  South  Jetties)  to 
the  estuary  and  the  Astoria-Megler  Bridge. 
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Figure  9. 


The  Mouth  of  the  Columbia  River  with  jetties,  Columbia  River  bar, 
and  Astoria-Megler  Bridge  labeled  (after  Google  Earth,  2014). 


This  region  was  chosen  because  this  is  the  portion  of  the  estuary  where  the  rivers 
strong  discharge  begins  to  enter  into  a  more  narrow  and  constricted  region  of  the  estuary. 
The  jetties  funnel  the  flow  through  in  a  jet-like  current  that,  during  ebb  currents,  can 
routinely  reach,  and  exceed,  2  m/s  (Gonzales  1984).  The  bar  (Figure  10)  that  has  formed 
outside  of  the  MCR  as  a  result  of  these  intense  wave-current  interactions  and  high 
volume  sediment  transport,  has  been  one  of  the  main  causes  of  approximately  2,000  ships 
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being  lost  at  sea  in  this  region  since  1792  (Simenstad  et  al.  2002).  The  bar  extends 
offshore  from  the  mouth  of  the  river  approximately  8  km.  Its  flat,  shallow  nature  in 
contrast  to  the  surrounding  waters  makes  navigating  this  bar  a  treacherous  ordeal.  The 
Navigation  Channel  has  been  dredged  through  the  southern  point  of  the  bar  leading  into 
the  MCR.  The  channel  depth  averages  17-20m  through  the  bar  and  continues  on  into  the 
MCR  to  bifurcate  into  the  North  and  South  Channels.  The  South  Channel  is  the  main 
navigation  channel.  At  approximately  -123.95°  longitude  there  is  a  trench  within  the 
South  Channel  that  is  approximately  20-25  m  deep  (almost  10  m  deeper  than  the 
surrounding  water). 
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Figure  10.  USGS  bathymetry  of  the  bar  and  the  MCR.  The  trench  is  well 
defined  in  the  main  navigation  channel  (after  Stevens  2014). 


19 


C.  SHIPBOARD  ACOUSTIC  DOPPLER  CURRENT  PROFILER 
ARRANGEMENT 


I,  Measurement  Principle 


The  Acoustic  Doppler  Current  Profiler  (ADCP)  is  the  most  commonly  used 
instrument  by  oceanographers  to  measure  currents  within  a  water  column.  The 
instrument  is  popular  for  a  multitude  of  reasons  to  include  a  wide  variety  of  operating 
frequencies  and  depths,  a  well-tested  and  documented  history  of  accuracy  and 
performance,  and  the  ease  with  which  the  instrument  can  be  deployed  in  various 
configurations  (i.e.,  hull-mounted,  over-the-side,  and  bottom-moored)  (RD  Instruments 
1996). 

ADCPs  work  by  using  the  Doppler  effect  to  measure  the  relative  radial  velocity 
between  different  objects  (Wewetzer  et  al.  1999).  The  Doppler  effect  is  the  observed 
change  in  frequency  of  a  wave  to  an  observer  moving  relative  to  its  source.  The  Doppler 
principle  can  be  described  using  a  water  wave  in  the  following  analogy.  A  stationary 
observer  (Figure  11a.)  is  sitting  on  the  surface  of  the  water  and  notices  that  waves  are 
passing  his  position  a  rate  of  one  wave  per  second.  If  the  observer  begins  to  move  in  the 
direction  of  the  source  of  the  waves  at  a  rate  of  four  waves  per  second  (Figure  I  Ib),  he 
will  perceive  a  passage  of  five  waves  per  second,  even  though  the  frequency  of  wave 
emission  has  not  changed  (Simpson  2001). 


Wave  rate: 


Wave  rate: 


Figure  11.  a.  Stationary  Observer 


b.  Moving  Observer 


(from  Simpson  2001) 
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ADCPs  utilize  the  Doppler  effect  by  transmitting  sound  at  a  fixed  frequency 
through  the  water  column  and  listening  for  the  echoes  returned  from  various  ‘scatterers’ 
or  suspended  particles  within  the  column  (RD  Instruments  1996).  The  particles  are 
assumed  to  be  of  such  size  and  strength  that  they  move  with  the  water  so  that  their 
velocities  are  representative  for  the  fluid  velocity. 

2,  Shipboard  ADCP  Configuration 

The  vertical  structure  of  ocean  currents  is  often  measured  with  an  ADCP  mounted 
to  the  hull  of  a  ship.  In  this  configuration  the  sensor  is  used  to  measure  the  movement  of 
particles  suspended  within  the  water  column  relative  to  the  ADCPs  position.  The 
navigation  system  is  used  to  remove  ship’s  motion  from  the  ADCP  data  during  post¬ 
processing  to  estimate  the  ocean  currents  in  a  fixed  reference  frame. 

The  R/V  Point  Sur  is  equipped  with  a  Teledyne  RD  Instruments’  Workhorse 
Mariner  ADCP  that  operates  at  300  kHz  with  a  typical  ping  rate  of  2  Hz  (Teledyne  2014). 
Depth  cell  size  for  this  data  was  set  at  2  m,  which  yields  a  range  of  78-102  m  with  a 
single  ping  standard  deviation  of  6.1  cm/s  (Teledyne  2014).  Further  information  on  the 
300  kHz  unit’s  performance,  accuracy,  and  resolution  can  be  found  at  the  following  link: 

http://marineops.mlml.calstate.edu/sites/default/files/wh  mariner  ds  Ir.pdf 

The  WH  300  ADCP  unit  is  mounted  to  the  bottom  of  the  R/V  Point  Sur  hull  at 
mid-ship.  The  unit  is  mounted  such,  that  it  can  be  accessed  both  internally  (transducer 
electronics  access)  and  externally  (biofouling  removal).  The  unit  is  recessed  so  as  to  not 
protrude  past  the  hull.  This  protects  the  ADCP  from  debris  damage,  and  minimizes  the 
effect  of  air  bubble  attenuation  across  the  path  of  the  ADCP  beams.  Attenuation 
(weakening)  of  the  signal  by  air  bubbles  reduces  the  profiling  range  and  thus 
performance  (RD  Instruments  2002). 

The  WH  300  ADCP  unit  utilizes  four  sound  beams  to  measure  current  velocity  in 
three  dimensions.  Data  is  divided  in  evenly  spaced  segments  throughout  the  water 
column  called  bins.  This  aspect  of  data  collection  makes  the  ADCP  very  similar  to  a 
moored  line  of  standard  current  meters.  The  ADCP  uses  its  four  transducers  to  measure 
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horizontal  velocity  as  well  as  vertical  velocity  of  the  currents  (Figure  12).  When  the 
instrument  beams  are  aligned  to  Earth  coordinates,  one  pair  of  beams  measures  the  north- 
south  component  of  velocity  while  the  other  pair  of  beams  measures  the  east-west 
component  of  velocity.  While  each  pair  of  beams  measures  a  different  component  of 
horizontal  velocity,  both  pairs  of  beams  are  simultaneously  measuring  vertical  velocity  as 
well  (RD  Instruments  1996).  This  is  where  ADCPs  differ  greatly  from  standard  current 
meters.  Taking  the  sum  of  opposing  beams  cancels  the  horizontal  velocities  leaving  only 
vertical  velocities.  Since  these  vertical  velocities  are  assumed  to  be  negligible  they  can 
thus  be  used  as  an  error  check  (Wewetzer  et  al.  1999). 


-  Beam  velocity'  component 


First  pair  of  beams 
calculates  east-west 
and  verical  velocity 


Second  pair  of  beams 
calculates  north-south 
and  vertical  velocity 


Figure  12.  The  relationship  of  beam  and  Earth  velocity  components 

(from  RD  Instruments  1996). 


3,  Post-Processing  Ship  Mounted  Data 

Data  collected  from  a  ship-mounted  ADCP  requires  the  correction  of  two  types  of 
motion:  rotation  (pitch,  roll,  and  heading),  and  translation  (ship  velocity).  Rotation 
correction  is  needed  because  the  ship’s  pitch  and  roll  affect  the  orientation  of  the  beams 
and  thus  the  velocity  components  that  are  being  measured.  These  distortions  are 
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compounded  with  an  increase  in  severity  of  surface  conditions  (RD  Instruments  1996). 
Translation  correction  is  also  essential  because  the  ship’s  velocity,  if  moving,  is  usually 
greater  than  the  actual  current  velocity  (RD  Instruments  1996). 

The  University  of  Hawaii  has  developed  a  suite  of  post-processing  codes  called 
UHDAS  (University  of  Hawaii  Data  Acquisition  System)  that  acquire  data  from  the 
ADCP  and  ancillary  sensors  such  as  GPS,  gyrocompass,  and  inertial  attitude  sensors  and 
use  CODAS  (Common  Ocean  Data  Access  System)  to  build  a  dataset  via  user  specified 
configurations.  CODAS,  through  a  series  of  editing  steps,  initially  uses  a  reliable 
heading  source  such  as  a  gyro  to  correct  for  heading  errors,  find  UTC  times  of  all  data 
points,  and  edit  out  bad  single-ping  velocities.  Finally,  the  reference  layer  is  found  and 
smoothed,  bottom-track/water-track  analysis  is  used  to  determine  preliminary  angle  and 
amplitude  corrections,  bad  data  caused  by  bottom  interference,  or  bubbles,  are  edited  out, 
notes  can  be  produced  for  later  use,  and  the  data  files  are  extracted  into  a  NetCDF  format 
for  use  in  MATLAB.  Further  information  on  the  UHDAS  Suite  and  CODAS  post¬ 
processing  method  can  be  found  at  the  following  website: 

http://currents.soest.hawaii.edU/docs/doc/index.html# 
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IV.  DATA  AND  MODEL  ANALYSIS 


A,  ADCP  TRANSECTS 

1,  Selection  Process 

To  investigate  the  vertical  structure  of  tidal  currents,  the  shipboard  ADCP  data 
was  broken  up  into  transects  where  the  ship  was  traversing  the  area  of  interest  at  a  nearly 
constant  heading.  The  complete  ship  track  shown  in  Figure  13  contains  numerous 
transects  in  the  shipping  channels  and  across  the  bar  when  the  ship  was  conducting  drifter 
deployments. 


20 1 3  MCR  Research  Cruise 


LongKudefdeg] 

Figure  13.  Complete  cruise  track  of  the  2013  RIVET  II  Research  Cruise. 
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The  first  step  involved  breaking  the  entire  data  set  down  into  transeets.  This  was 
aehieved  by  an  automated  termination  of  transects  when  the  ship  performed  a  change  in 
heading  greater  than  a  user-specified  maximum.  The  ADCP  logged  profiles  at  two- 
minute  intervals,  and  heading  changes  were  measured  between  successive  profiles.  The 
transects  were  terminated  when  the  heading  changed  more  than  30  degrees.  To  consider 
only  transects  with  sufficient  duration  in  the  analysis,  transects  shorter  than  40  minutes 
were  discarded. 

Next,  transects  were  manually  checked  for  their  geographical  location.  As  seen  in 
Figure  13,  most  of  the  data  collected  was  collected  outside  of  the  MCR.  Only  transects 
that  were  located  outside  of  the  river  in  close  proximity  to  the  bar,  and  inside  the  MCR 
between  the  Astoria-Megler  Bridge  and  the  jetties,  were  kept.  Any  transects  located 
outside  of  the  MCR  away  from  the  bar,  and  too  far  upriver  were  discarded. 

Finally,  all  transects  were  manually  checked  for  speed.  Large  decreases  in  ship’s 
speed  could  possibly  indicate  the  ship  holding  station.  As  with  heading,  speed  was 
calculated  using  successive  points  within  the  data.  Transects  that  exhibited  changes  in 
speed  greater  than  3  m/s  or  near  zero  speed  were  discarded. 

2,  Selected  Sections 

Once  the  selection  process  was  complete,  vertical  section  plots  were  generated 
that  display  the  current  magnitude  and  direction  with  longitude  vs.  depth.  Of  the  16 
transects  (Figure  14)  that  were  kept,  11  of  them  were  located  in  the  main  shipping 
channel  between  the  mouth  of  the  river  and  the  Astoria-Megler  Bridge.  The  remaining 
five  transects  are  offshore  over  the  bar  area  and  its  approaches 

Predicted  tidal  currents  were  retrieved  from  NOAA  archives  for  the  Clatsop  Spit 
station.  This  data  contains  the  predicted  magnitude  and  times  for  peak  ebb  and  flood 
currents  as  well  as  slack  water.  The  start  and  ending  times  for  each  transect  were  used  in 
conjunction  with  the  predicted  tidal  current  data  to  identify  what  phase  of  the  tidal  cycle 
each  transect  was  occurring  in. 
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Figure  14.  Image  displaying  the  locations  of  all  transects  (white  lines) 

(after  Google  Earth  2014). 


Transects  were  initially  divided  according  to  their  geographical  location;  inside 
the  MCR  and  outside  the  MCR.  Next,  all  transects  were  further  divided  into  transects 
occurring  during  flood,  and  transects  occurring  during  ebb.  Finally,  transects  were 
further  divided  according  to  their  proximity  to  slack  water: 

•  Transects  that  have  experienced  slack  water  and  are  approaching  peak 
flood. 

•  Transects  that  have  experienced  peak  flood  and  are  approaching  slack 
water. 

•  Transects  that  have  experienced  slack  water  and  are  approaching  peak 
ebb. 

•  Transects  that  have  experienced  peak  ebb  and  are  approaching  slack 
water. 

Transects  were  divided  in  this  manner  because  one  of  the  goals  of  this  research 
was  to  characterize  the  vertical  structure  of  the  tidal  currents  during  different  tidal  stages. 

B,  SELFE  MODEL  DATA 

To  help  limit  the  size  of  the  model  data  files,  SELFE  model  data  was  extracted 
based  on  the  times  and  locations  of  the  observational  transects.  Once  all  the  model  data 


27 


were  extracted  for  each  ADCP  transect,  a  model  transect  was  constructed  by  choosing  the 
nearest  model  grid  point  in  space  and  time  for  each  ADCP  profile. 


C.  FINAL  SELECTION 

Comparisons  between  the  observational  data  and  the  model  were  made  for  all  16 
transects.  A  subset  of  these  transects  that  characterized  the  vertical  structure  at  each  tidal 
stage  are  presented  in  the  next  two  chapters.  Transects  were  purposely  chosen  during 
early  and  late  flood  and  ebb  periods  to  illustrate  the  difference  in  flow  structure  in 
different  parts  of  the  tidal  cycle.  These  final  six  transects  (four  inside  the  river  and  2 
outside  the  river)  are  shown  in  Figure  15. 
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Figure  15.  Final  six  transects,  (a)  Shows  all  six  transects,  (b)  Shows  four 
transects  within  the  MCR.  (c)  Shows  two  transects  outside  the  MCR 
close  to  the  bar  (after  Google  Earth  2014). 
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V.  FLOOD  RESULTS 


Data  collection  for  the  ADCP  begins  at  a  depth  below  the  surface  of  7.08  m  due 
to  ADCP  blank  range.  This  blank  range  serves  multiple  purposes  by  being  a  minimum 
distance  necessary  to  avoid  collecting  bad  data  that  has  potentially  been  contaminated 
from  ringing,  and  avoid  flow  noise  from  the  bottom  of  the  ship  when  hull-mounted  (RD 
Instruments  2002). 

Each  transect  observation  and  model  vertical  section  are  displayed  for  comparison 
along  with  transect  start  and  end  times,  and  NOAA  predicted  tidal  current  magnitudes 
and  times.  Each  vertical  section  displays  depth,  longitude,  and  current  speed  and 
direction.  Additionally,  current  vectors  are  plotted  for  selected  depths  to  allow  a  more 
quantitative  comparison  of  observation  and  model  data.  Transects  will  be  divided  for 
analysis  by  flood  (Eigure  16)  and  ebb  (Figure  22). 
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Figure  16.  Transects  during  flood  inside  the  MCR  (left  panel)  and  outside  the 

MCR  (right  panel). 
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A.  TRANSECT  1:  PEAK  FLOOD  TO  SLACK  INSIDE  MCR 

Transect  1  was  recorded  on  06  Jun  2013  between  1946-2046  UTC.  The  ship  was 
heading  southeast  upriver  and  peak  flood  was  predicted  to  occur  at  1838  UTC  with  slack 
predicted  to  occur  at  2114  UTC.  New  moon  occurred  on  08  Jun  just  two  days  after, 
indicating  this  transect  was  very  close  to  a  neap  to  spring  transition. 

Observational  data  (Figure  17)  shows  that  the  current  is  directed  upriver  from  the 
uppermost  bin  (7.08  m  below  the  surface)  throughout  the  remainder  of  the  column.  The 
effects  of  the  bottom  boundary  layer  are  evident  in  the  current  weakening  with  depth 
towards  the  bottom.  Stronger  currents  at  the  beginning  of  the  period  drastically 
weakened  from  86  cm/s  to  approximately  20  cm/s  in  the  middle  of  the  period.  Tidal 
current  strength  unexpectedly  increased  toward  the  end  of  approaching  slack,  predicted  to 
occur  just  28  minutes  after  the  end  of  the  transect. 

SELFE  model  data  shows  a  somewhat  different  picture.  The  strongest  currents 
were  predicted  to  be  at  the  start  of  the  period,  as  with  the  observational  data,  and  the 
observed  and  predicted  initial  current  speeds  are  in  good  agreement.  As  the  period 
progressed  toward  slack  water  the  model  shows  currents  weakening  throughout  the 
column  and  even  changing  direction  at  the  surface  before  the  end  of  the  period,  indicating 
the  development  of  an  ebbing  current  at  the  surface  while  the  lower  water  column  is  still 
flooding.  While  this  phenomenon  is  well  established,  it  could  not  be  verified  with  the 
ADCP  data,  which  do  not  cover  the  upper  seven  meters  of  the  water  column.  Towards 
the  end  of  the  transect  the  model  predicts  notably  weaker  currents  near  the  bottom  than 
observational  data. 

A  closer  comparison  of  this  data  using  current  vector  plots  (Eigure  18)  at  7.08  m 
and  13.08  m  for  both  sets  of  data  shows  differences  as  well.  At  7.08  m,  the  observed 
current  at  the  beginning  of  the  period  has  a  strong  eastern  component  before  gradually 
turning  south  and  heading  upriver.  The  current  maintains  this  southern  flow  until  - 
123.94°  longitude;  in  the  vicinity  of  the  trench  identified  in  Eigure  10.  It  is  possible  that 
this  trench  is  topographically  steering  the  current  back  to  an  eastern  direction  of  flow. 
Model  data  at  7.08  m  shows  a  predominantly  southern  current  along  the  entire  transect. 
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The  model  data  shows  the  eurrent  weaken  much  more  toward  the  end  of  the  period  than 
does  the  observational  data. 


Transect  Times 


06Jun  1946 
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Figure  17.  Transect  1  vertical  section  plots  for  ADCP  (top  panel)  and  SELFE 
(bottom  panel)  data.  Current  speed  is  shown  in  color  with  positive 
values  indicating  up-river  (flood)  flow  and  negative  values  indicating 
down-river  (ebb)  flow.  NOAA  predictions  of  max  flood  current  and 
times  of  max  flood  and  slack  water  at  Clatsop  Spit  station  are  indicated 

on  the  left. 
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ADCP  Current  Vectors  at  7.08m  depth 


ADCP  Current  Vectors  at  13.08m  depth 


Longitude 


Longitude 


SELFE  Current  Vectors  at  7.08m  depth 


Longitude 


SELFE  Current  Vectors  at  13.08m  depth 


Longitude 


Figure  18.  Transect  1  current  vectors  for  ADCP  (top  panels)  and  SELFE 

(bottom  panels)  data  at  7.08  m  (left  panels)  and  13.08  m  (right  panels) 

depth. 


B.  TRANSECT  2:  SLACK  TO  PEAK  FLOOD  INSIDE  MCR 

Transect  2  was  recorded  on  02  Jun  2013  between  1340-1424  UTC.  The  ship  was 
heading  northwest  out  to  sea  and  slack  was  predicted  to  occur  at  1145  UTC  with  peak 
flood  predicted  to  occur  at  1508  UTC.  Third  quarter  moon  occurred  on  31  May, 
indicating  this  transect  was  recorded  during  a  neap  cycle. 

ADCP  vertical  section  plots  (Figure  19)  show  weaker  currents  at  the  beginning  of 
the  transect  after  the  passage  of  predicted  slack  water.  The  strongest  currents  occur  in  a 
region  in  the  center  of  the  transect  that  eventually  begin  to  weaken  as  the  transect  period 
approaches  predicted  peak  flood,  which  was  predicted  to  occur  44  minutes  after  the  end 
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of  the  transect.  Again,  this  was  an  unexpected  result.  Assuming  the  predicted  tidal 
current  times  to  be  accurate,  it  is  reasonable  to  assume  the  tidal  currents  should  begin  to 
increase,  especially  considering  how  close  the  end  of  the  transect  is  to  peak  flood. 
Instead  the  observation  data  shows  peak  flood  occurring  in  the  middle  of  the  transect 
period  and  by  the  end  is  already  transitioning  to  slack  again.  At  14  m  the  tidal  currents 
begin  to  approach  zero  and  even  shift  direction  to  flow  out  of  the  river.  Prior  to  this, 
however,  the  currents  strength  is  increasing  with  depth. 

SELFE  predictions  agree  with  the  ADCP  observed  weaker  currents  at  the 
beginning  of  the  transect  period  and  the  region  of  strong  current  in  the  center  of  the 
transect.  Model  data,  however,  does  show  stronger  currents  toward  the  end  of  the  period 
as  the  transect  approaches  predicted  peak  flood.  SEEFE’s  use  of  a  tide  model  that 
depends  on  predicted  tides  itself  could  be  a  possible  explanation  for  this  difference.  The 
model  data  and  observational  data  also  disagree  at  the  bottom  boundary.  The  model  does 
not  show  the  current  approaching  zero  as  quickly  as  the  observational  data  does.  The 
model  data  also  does  not  show  current  reversal  at  the  bottom. 

Further  examinations  into  this  transect  using  current  vector  plots  (Figure  20) 
reveals  good  agreement  among  the  observations  and  model  with  regards  to  direction  and 
magnitude  of  the  current  throughout  the  water  column.  At  13.08  m  both  observations  and 
model  show  a  sudden  southward  shift  in  the  current  direction  in  the  vicinity  of  the 
previously  discussed  trench.  Again,  this  appears  to  be  topographic  steering  occurring 
because  of  the  trench,  but  more  importantly,  the  model  clearly  predicts  this  behavior. 
The  observational  data  shows  a  clear  reversal  at  the  end  of  the  transect  period  at  13.08  m 
that  is  not  predicted  by  the  model. 
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Figure  19.  Transect  2  vertical  section  plots  for  ADCP  (top  panel)  and  SELFE 
(bottom  panel)  data.  Current  speed  is  shown  in  color  with  positive 
values  indicating  up-river  (flood)  flow  and  negative  values  indicating 
down-river  (ebb)  flow.  NOAA  predictions  of  max  flood  current  and 
times  of  max  flood  and  slack  water  at  Clatsop  Spit  station  are  indicated 

on  the  left. 
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ADCP  Current  Vectors  at  7.08m  depth 


Longitude 


ADCP  Current  Vectors  at  13.08m  depth 


Longitude 


SELFE  Current  Vectors  at  7.08m  depth 


Longitude 


SELFE  Current  Vectors  at  13.08m  depth 


Longitude 


Figure  20.  Transect  2  current  vectors  for  ADCP  (top  panels)  and  SELFE 

(bottom  panels)  data  at  7.08  m  (left  panels)  and  13.08  m  (right  panels) 

depth. 


C.  TRANSECTS:  FLOOD  OUTSIDE  OF  MCR 

Transect  5  was  recorded  on  03  Jun  2013  between  0124-0214  UTC.  The  ship  was 
heading  southeast  toward  the  river  mouth  and  slack  was  predicted  to  occur  at  02  Jun  2013 
2315  UTC  with  peak  flood  predicted  to  occur  at  0314  UTC.  Third  quarter  moon 
occurred  on  3 1  May  indicating  this  transect  was  recorded  during  a  neap  cycle. 

Vertical  sections  (Figure  21)  for  this  transect  were  produced  for  both  the  u  and  v 
components  of  velocity.  Vertical  sections  for  the  u  component  of  velocity  show  a  mild 
shoreward  current  throughout  the  column  with  weak  seaward  flow  near  the  bottom. 
Vertical  sections  for  the  v  component  of  velocity  show  a  weak  southern  current 
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throughout  the  region  with  some  intensifieation  near  the  bottom  in  the  middle  of  the 
transeet. 

Model  data  for  both  eomponents  of  veloeity  shows  homogeneous  eurrents 
throughout  the  water  eolumn,  with  the  exception  of  the  upper  few  meters,  and  do  not 
depict  any  the  vertical  structure  with  flow  reversals  that  observations  reveal.  The  models 
inability  to  define  structure  outside  of  the  MCR  is  because  the  MCR  SELFE  model  grid 
is  tuned  to  focus  on  the  estuary.  The  model  grid  outside  of  the  estuary  is  coarser,  which 
negatively  affects  accuracy.  This  region  also  lacks  the  substantial  observational  network 
that  the  estuary  possesses  (Kama  2014). 


03Jun  0124 


03Jun  0214 


02Jun  2315  sufj  u  v*tocitv  (cm/s| 


l4n0lwde 


Figure  21.  Transect  5  vertical  section  plots  for  ADCP  observations  (upper 
panels)  and  SELFE  predictions  (lower  panels).  Left  panels:  u 
component  (+  eastward,  -  westward).  Right  panels:  v  component 
(+  northward,  -  southward) 
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VI.  EBB  RESULTS 
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Figure  22.  Transects  on  the  left  represent  ebb  data  inside  the  MCR  while  the 
transect  on  the  right  represents  ebb  data  outside  of  the  MCR. 


A,  TRANSECTS:  PEAK  EBB  TO  SLACK  INSIDE  THE  MCR 

Transect  3  was  recorded  on  09  Jun  2013  between  0240-0320  UTC.  The  ship  was 
heading  southeast  upriver  and  peak  ebb  was  predicted  to  occur  at  0050  UTC  with  slack 
predicted  to  occur  at  0351  UTC.  New  moon  occurred  on  08  Jun,  indicating  this  transect 
was  recorded  during  a  spring  cycle. 

Slack  water  was  predicted  to  occur  31  minutes  after  the  end  of  the  transect, 
preceding  the  transition  back  to  flood.  ADCP  vertical  section  (Figure  23)  plots  show  the 
flood  current  already  forming  at  the  beginning  of  the  transect  period.  Strong  ebb  currents 
remain  in  the  upper  part  of  the  water  column  and  continue  to  ebb  throughout  the  transect 
period.  At  10  m  the  current  strength  approaches  zero  and  begins  to  flood  around  13  m. 
This  is  an  interesting  result  because  not  only  is  the  flood  transition  occurring  well  before 
the  predicted  times,  but  the  strong  flood  currents  near  the  bottom  while  the  surface  waters 
are  still  ebbing  indicates  there  is  no  slack  at  all. 
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SELFE  model  predictions  agree  with  the  overall  structure  and  timing  of  these 
events,  but  with  some  disparity  in  magnitude.  Model  data  also  shows  the  current 
continuing  to  ebb  over  the  entire  water  column  with  weak  flows,  indicating  a  boundary 
layer  near  the  bottom.  The  model  also  shows  a  period  of  time  at  the  end  of  the  transect  in 
which  the  entire  column  approaches  slack  water.  Weak  upriver  currents  are  predicted  at 
the  bottom,  but  overall  the  model  is  showing  a  transition  to  slack  water. 

Further  analysis  using  current  vectors  (Figure  24)  shows  good  agreement  between 
the  observations  and  the  model  at  7.08  m.  Both  observations  and  model  prediction  data 
show  an  ebbing  current  throughout  the  period  at  this  depth,  although  the  observational 
data  shows  a  stronger  current.  Observations  and  model  predictions  differ  greatly  at  13.08 
m.  Observational  data  shows  the  current  transitioning  into  a  strong  flood  current  at 
approximately  a  third  of  the  distance  of  the  transect.  Model  data  shows  the  current 
weakening  and  approaching  zero,  and  transitioning  into  a  much  weaker  flood  current 
towards  the  end  of  the  period.  Again,  the  model  current  vectors  at  this  depth  resemble  a 
transition  to  slack  water,  whereas  the  observational  data  at  this  depth  clearly  show  a 
transition  to  flood  current. 
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Figure  23.  Transect  3  vertical  section  plots  for  ADCP  (top  panel)  and  SELFE 
(bottom  panel)  data.  Current  speed  is  shown  in  color  with  positive 
values  indicating  up-river  (flood)  flow  and  negative  values  indicating 
down-river  (ebb)  flow.  NOAA  predictions  of  max  ebb  current  and 
times  of  max  ebb  and  slack  water  at  Clatsop  Spit  station  are  indicated 

on  the  left. 
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ADCP  Current  Vectors  at  7.08m  depth 


ADCP  Current  Vectors  at  13.08m  depth 
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SELFE  Current  Vectors  at  7.08m  depth 
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Figure  24.  Transect  2  current  vectors  for  ADCP  (top  panels)  and  SELFE 

(bottom  panels)  data  at  7.08  m  (left  panels)  and  13.08  m  (right  panels) 

depth. 


B.  TRANSECT  4:  SLACK  TO  PEAK  EBB  INSIDE  THE  MCR 

Transect  4  was  recorded  on  30  May  2013  between  1458-1558  UTC.  The  ship 
was  heading  southeast  upriver  and  slack  was  predicted  to  occur  at  1344  UTC  with  peak 
ebb  predicted  to  occur  at  1644  UTC.  Third  quarter  moon  occurred  on  31  May,  just  one 
day  after  the  recording  of  this  data,  indicating  this  transect  was  transitioning  from  spring 
to  neap  cycle. 

ADCP  vertical  sections  (Eigure  25)  show  a  strong  (2  m/s)  ebb  current,  with  a 
maximum  at  the  upper  most  bin  (7.08  m)  and  much  weaker  currents  near  the  bottom.  As 
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the  transect  progresses  the  stronger  currents  extend  further  down  the  column,  and  by  the 
end  the  period  a  strong  well-defined  ebb  current  is  observed  at  all  depths. 

SELFE  model  predictions  are  in  excellent  agreement  with  the  observational  data. 
There  is  a  strong  ebb  current  throughout  the  column  with  weaker  currents  in  the  center  of 
the  transect  near  the  bottom.  The  region  of  weakening  current  is  smaller  in  the  model 
data,  but  is  represented  in  the  same  location.  Observed  and  predicted  current  strength 
agree  well. 

Further  analysis  using  current  vectors  (Figure  26)  shows  good  agreement  between 
observations  and  model  predictions  At  13.08  m  depth,  observed  and  predicted  current 
direction  show  the  current  slightly  veering  to  the  north  around  the  vicinity  of  the  trench 
identified  within  the  main  navigation  channel.  This  shift  to  north  vice  south  during  flood 
is  consistent  with  a  topographic  steering  effect.  Overall,  observed  and  predicted  vertical 
current  structures  are  in  excellent  agreement  on  this  transect. 
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Figure  25.  Transect  4  vertical  section  plots  for  ADCP  (top  panel)  and  SELFE 
(bottom  panel)  data.  Current  speed  is  shown  in  color  with  positive 
values  indicating  up-river  (flood)  flow  and  negative  values  indicating 
down-river  (ebb)  flow.  NOAA  predictions  of  max  ebb  current  and 
times  of  max  ebb  and  slack  water  at  Clatsop  Spit  station  are  indicated 

on  the  left. 


44 


ADCP  Current  Vectors  at  7.08m  depth 
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Figure  26.  Transect  4  current  vectors  for  ADCP  (top  panels)  and  SELFE 

(bottom  panels)  data  at  7.08  m  (left  panels)  and  13.08  m  (right  panels) 

depth. 


C.  TRANSECT  6:  EBB  OUTSIDE  OF  MCR 

Transect  6  was  recorded  on  31  May  2013  between  1900-2048  UTC.  The  ship 
was  heading  southeast  toward  the  river  mouth  and  peak  ebb  was  predicted  to  occur  at 
1750  UTC  with  slack  predicted  to  occur  at  2121  UTC.  Third  quarter  moon  occurred  on 
3 1  May  indicating  this  transect  was  recorded  during  a  neap  cycle. 

Vertical  sections  (Figure  28)  for  this  transect  were  produced  for  both  the  u  and  v 
components  of  velocity.  Vertical  sections  for  the  u  component  of  velocity  show  a  weak 
seaward  current  on  the  surface  with  a  stronger  shoreward  current  in  the  bottom  20  m. 
Vertical  sections  for  the  v  component  of  velocity  show  a  strong  northward  current  in  the 
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top  20  m  the  plot.  Below  this  depth  the  current  weakens  and  turns  to  flow  southward 
along  the  bottom. 

Model  predictions  show  relatively  homogeneous  conditions  throughout  the 
column  and  stronger  current  in  the  upper  10  m  of  the  water  column.  The  model  does 
show  a  strong  ebb  current  near  the  surface  in  the  second  half  of  the  transect  (negative  u, 
positive  v)  in  qualitative  agreement  with  observations  but  the  detailed  vertical  structure 
looks  rather  different. 
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Figure  27.  Transect  6  vertical  section  plots  for  ADCP  observations  (upper 
panels)  and  SELFE  predictions  (lower  panels).  Eeft  panels:  u 
component  (+  eastward,  -  westward).  Right  panels:  v  component 
(+  northward,  -  southward) 
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VII.  CONCLUSION 


A.  SUMMARY 

1.  MCR  Shipboard  ADCP  Measurements 

The  vertical  tide  structure  of  the  MCR  is  not  well  understood,  in  large  part,  due  to 
a  lack  of  sufficient  measurements.  The  objectives  of  this  research  were  to  characterize 
the  vertical  structure  of  the  tidal  within  the  MCR  over  various  tidal  cycles,  and  validate 
the  SELFE  model’s  ability  to  predict  these  currents. 

To  investigate  the  vertical  structure  of  the  tidal  currents,  ADCP  measurements 
collected  during  a  three-week  long  cruise  within  the  MCR  were  broken  up  into  transects. 
This  was  accomplished  by  extracting  periods  of  data  collected  while  traversing  the  main 
navigation  channel  and  the  bar  region  at  a  relatively  constant  speed  and  heading.  The 
initial  three-week  data  set  was  reduced  to  16  transects:  11  inside  the  river  and  five 
outside  the  river,  seaward  of  the  river  mouth.  A  subset  of  these  transects  are  presented 
here  that  characterize  the  vertical  structure  at  different  tidal  stages. 

2,  Observed  Vertical  Structure 

Early  in  the  flood  cycle  the  vertical  structure  exhibits  two-layer  flow  with  a  strong 
flood  current  along  the  bottom  and  a  slight  river  outflow  on  the  surface.  As  the  period 
approaches  peak  flood,  the  entire  water  column  experiences  a  complete  reversal  in  flow 
to  a  flood  current.  This  complete  reversal  in  flow  through  the  entire  water  column  can  be 
seen  as  far  upriver  as  the  Astoria-Megler  Bridge. 

The  ebb  current  is  more  homogenous  throughout  the  water  column  than  the  flood 
current.  River  outflow  acts  in  unison  with  the  ebb  current  to  produce  a  stronger  current. 
During  higher  range  tidal  periods  slack  water  is  often  bypassed  altogether,  with  a  flood 
current  developing  in  the  lower  part  of  the  water  column  while  strong  ebb  flows  continue 
in  the  surface  layer.  This  sudden  transition  from  one  strong  tidal  cycle  to  the  next  with 
no  slack  water  in  between  is  one  of  the  many  dangers  mariners  face  when  navigating  the 
channel. 
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3. 


SELFE  Model  Evaluation 


The  SELFE  model  is  an  operational  flow  model  developed  and  run  by  the  Oregon 
Health  and  Seience  Elniversity.  The  SELFE  model  solves  the  shallow-water  equations  on 
an  unstruetured  grid,  whieh  allows  the  resolution  to  vary  within  a  single  model  domain. 
The  model  solves  its  primary  variables  in  three  dimensions.  Model  inputs  eome  from 
observation  stations  within  the  river,  meteorological  databases,  and  oceanographic 
forcing  models. 

Inside  the  river,  the  SELFE  model  performed  qualitatively  well  in  comparison  to 
the  observations.  The  model  showed  good  agreement  for  nine  out  of  the  original  11 
transects  inside  the  river,  whereas  Transects  1  and  3,  showed  significant  differences 
between  model  and  observations.  The  data  from  these  two  transects  were  collected 
during  periods  approaching  slack.  In  the  transition  from  flood  to  slack  the  predictions 
show  weaker  currents  in  the  lower  part  of  the  water  column  than  is  observed,  whereas  in 
the  transition  from  ebb  to  slack  the  predicted  reversal  to  a  near-bottom  flood  current  is 
weaker  than  is  observed. 

Outside  of  the  river  the  SELFE  model  does  not  resolve  the  detailed  vertical 
structure  of  the  observed  currents.  This  is  not  surprising  as  the  grid  resolution  outside  of 
the  mouth  is  coarser  than  the  grid  resolution  inside  the  mouth.  The  area  of  the  model 
domain  outside  of  the  mouth  of  the  river  also  suffers  from  limited  observational  inputs 
whereas  the  estuary  has  a  well-established  observational  network  to  which  the  model  is 
calibrated. 

B,  FUTURE  RESEARCH 

The  ADCP  is  becoming  a  standard  instrument  included  in  many  ships’ 
inventories.  Collecting  ADCP  data  while  vessels  transit  through  coastal  and  tidal  inlets 
could  help  shed  light  on  the  vertical  structure  and  wave-current  interactions  occurring  in 
these  dynamic  regions.  These  observations  are  useful  in  the  immediate  evaluation  of 
local  area  models  and  dynamics.  Collecting  these  observations  and  storing  them  in  a 
database  can  also  allow  them  to  be  used  in  the  future  for  characterizing  the  dynamics  in 
specific  regions  around  the  world. 
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The  development  of  an  automated  analysis  seheme  eould  make  the  collection  and 
evaluation  of  these  data  simpler  for  mariners  who  do  not  possess  the  knowledge  and 
expertise  required  to  manage  these  data.  This  would  allow  the  data  to  be  collected  by  a 
crew,  stored  during  their  cruise,  and  delivered  upon  return  to  port  to  a  central  location  for 
automated  analysis  and  archiving. 

The  ADCP  observations  collected  during  this  experiment  do  not  provide 
information  in  the  top  7  m  of  the  water  column,  thus  excluding  comparison  of  surface 
observations  with  the  model  surface  predictions.  Surface  drifter  data,  however,  were  also 
collected  in  the  experiment  (Pearman  2014).  Analysis  of  the  combined  ADCP  and  drifter 
data  will  allow  for  a  more  complete  analysis  of  the  vertical  structure  of  the  tidal  currents 
in  the  MCR. 
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